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CONVERGENT  WAVE  IN  A  PLASTIC 
MEDIUM 


/Following  is  the  translation  of  an  article 
by  E.  I.  Andriankln  entitled  Skhodyashch- 
ayasya  Volna  v  Plasticheskoy  Srede  (Engl¬ 
ish  version  above )  in  Doklady  Akademii  Nauk 
SSSR  (Proc.  Acad.  Sci.  USSR)  Vol.  No.  4 
If 4%  pages  659-662/  n  1 


(Presented  by  Academician  N.  N .  Semenov ,  13  Nov¬ 
ember  1959) 


The  problem  of  an  explosion  In  a  plastic  compres¬ 
sible  medium  has  been  Investigated  by  A.  S.  Kompaneyets 
in  /!/.  In  this  work  a  simplified  law  of  compression  was 
assumed.  This  made  It  possible  to  obtain  simple  equations 
for  the  law  of  motion  of  the  front  of  a  diverging 

spherical  wave.  , 

Utilizing  the  analogous  law  ox  compression  (we  will 
assume  that  the  density  of  the  medixxm  to  aches  the  limit 
yalue  p#  Pq  at  any  pressure  different  from  zero), 

we  will  consider  the  problem  of  the  converging  plastic 
wave,  which  can  be  formulated  as  follows. 

At  the  Instant  t  =  XQ,  let  there  be  a  pressure  at 

the  free  surface  of  a  .  spherical  layer  with  an  initial 
autside  radius  aQ.  Let  the  said  pressure  vary  according 

bo  a  given  pattern  P  *  PQF(t/X0)  (in  a  more  general  case 

i?  —  p(t/XQ,  x)  can  be  assumed).  At  the  Inner  surface  of 

the  spherical  layer  of  initial  radius  bQ  the  pressure  is 

zero.  At  the  instant  the  pressure  Is  applied,  a  shock 
•rave  is  propagated  through  the  material. 

We  assume  the  medium  behind  the  front  to.  be  incom- 


pressible,  •  p~~  and  the  condition  of  plasticity 

to  be  satisfied: 


~f*  Ifi  -4-  Qsq}, 


/  <  * 
0)  S 


where  k  and  m  are  assumed  to  be  known  constants,  dr  and 

6<p  are  the  main  stresses.  The  subscript  f  denotes 
the  values  at  the  wavefront.  Further,  we  assume 

| that  the  initial  density  of  the  material  in  the  spherical 
{layer  depends  on  the  radius  Pq"  Pq(vq)* 

We  shall  solve  the  problem  in  Lagrangian 

variables,  where  the  equations  of  continuity  will  be: 


dr  'oPo (ro) 


dr  dcr 


:  <ar  ~  °») 


Here  u  —  c)r/  c)t  -  r,  t  is  time,  r  and  r~  are  the  variable 
and  initial  coordinates  of  the  particles. 

Using  (1)  and  (2),  we  can  bring  Eq.  (3)  into 
the  form: 

ik[r' (i  -/J)]  =-■  ?»  ^o) rr~%  w  •  (4> 

'  p=~'" 

The  boundary  conditions  for  Eqs.  (2)  and  (4) 
are  the  laws  of  conservation  at  the  wavefront,  the  equal¬ 
ity  of  pressure  at  the  free  surface  rn  —  a  ,  and  the  con¬ 
dition  of  continuity  of  the  variable  Radius 

P*  ~  H  (R)  s^R2,  -  e+/Jt  r4>  =  R  «  r#f 

‘  (a° '  P  ~  P<jl'  (//to),  £*  =  1  —  po/p,-  ^ 

We  shall  introduce  the  dimensionless  parameters 
a  ~  aja„.  b  -  bja„,  x  —  R/a„.  s  =  r„/a0.  r  ~  rfu,t. 
o  —  p/o,.  xsrk/3mP0 ,  k  =  k/Rf,.  %  —  tj ?».  u  =Jrjdt,  {(>> 
dxfd~  ~  —Vy,  y  --  f >.R2/Pv  p.  ~  au  y 
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I 


_ _ _ _ _ _ _ I . — - - 

where  a*  and.  To*  are  the  outside  and  inside  radii  of  the 
spherical  layer  at  the  instant  t. 

Integrating  (2)  and  using  Conditions  (5b  we 
shall  find  s  i 

r;l  -=  .r'4~  3  \  (s)  s2</s,  a'1  =  ,v::  -f-  3  \  o„s-ds, 

*  _  *  _  (7) 

u  -  —  ,  x  =  z^x-  V y  4) .  s*  -  i  -  p..  4)- 


Integrating  Eq.  (4)  and  using  (7b  we  get 

xa)  -f  .y’P (*,  —  V7/  \  po  (h  r°-As-ds  -J-  2X2  \  p„  V  ?s-(/s, 


«=  x.  (r* 


bi>  =  Pi,  4)  a  =4(>.  4) 

If  the  law  of  motion  of  the  wavefront  is  known,  the 
pressure  and  velocity  distribution  in  the  entire  region 
x4  s  ^1, 

An  ordinary  differential  equation  for  the  vel¬ 
ocity  of  the  wavefront  can  he  obtained  from  (8),  if  we 
use  the  condition  on  the  free  surface.  Assuming  also 
C'  ~~  -f5,  sn,  while  calculating  all  integrals  in  Eqs,  (7) 
ana  (8).  we  shall  find: 


<h 

dx 

Here 


K  4)  z  4-  A’  4, '), 


dy 

dx 


I  -  P>A"\ 


xH(x) 


,  t  --  K 


*  =  1*  (9) 


P'  Pi,  4»)  “=■  Pp,, 


/c»— 


(a  -  -  -i i  (a*'  1  — xtt  1  ) 


)  ,  2(3  (a  —  I)  .v"*!-3-2 

-  *r  - 


r-+3ill=4il 

n  + 
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In  order  to  integrate  Eq,  (9)  it  is  necessary  to 
enow  its  asymptotic  expression  at  x  0  and  the  value  of 
the  derivative  dz/dx  at  point  x  -1* 


) 

i 

i 
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After  a  series  of  calculations*  we  shall. find. 


=*(£), 

-te, 

+  8- 

-  V (jr  K  +  73  U  *  {'  1)1 

(10) 

Unit 

x 

=  /i. 

t  »  -p . . .  (x-^Q,  n  >  0); 

(11) 

.2  =r  H  >x 

“  +  •  • 

•  ^ 

*<0). 

at  — ~  1 

(12) 

I  However*  if  the  outside  pressure  drops  rapidly*  the 

wave  may  stop  before  reaching  the  center. 

Equations  (8)  and  (9)  are  true  in  caseoC^Oj  if 
m  =0*  the  condition  of  plasticity  becomes  <<*,  -  <50  -  k  and 

the  solution  is  given  by  the  same  equations*  in  which  we 
have  to  assume 


*  ==  0,  oi  =b  2k  In  (x/a), 


It  should  also  be  noted  that  Eqs.  (9)  give  a  sol¬ 
ution  to  the  problem,  of  the  collapse  of  a  spherical 

layer  of  incompressible  plastic  material.  In  that  case  it 
is  necessary  to  assume 


r'  — -  s:!  -f-  x  -  b.  7\  =•  0, 

4(g-  1  ^  ia*  *  2  ( a  - 

(a -4)  (a*-  '  —  )  '*  A  ~ 


2  (a  —  1)  a-  fo,  — 

/.*  —  t  _  v=t—  1 


In  these  cases  an  asymptotic  solution  is  given  by 
Iq.  (11).  At  qr  -  1  and  x  -  bQ* 

the  boundary  conditions  for  this 
problem  will  be  y  ^  0;  in  the  case  of  initial  motion  in 
;he  fluid*  the  boundary  condition  will  be  y  -  The 
Latter  occurs*  for  instance*  when  a  spherical  layer  of 
jompressible  material  has  been  compressed  by  the  outside 
Dressure  to  the  limit  condition  p,,  —  f*:  then,  after 
bhe  disappearance  of  the  shock  wave  (at  the  boundary  with 
the  empty  space)  the  compressed  material  continues  to 
jollapse  as  if  Incompressible. 

In  Figs.  1  and  2  the  results  of  the  numerical  in¬ 
tegration  of  Eqs. (8)  are’  shown.  Curve  1  corresponds  to 
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the  solution  of  the  problem.,  for  parameter  values  of 


F  srt 


e**o.i 


n  -  0 


X  -  -o.i 


for  Curve  2.,  we  have 


P  r  t' 


£*-  0*1 


n  r  0 


k  r  -0.1 


fe  for  Curve  3>  we  have 


€  *  0,2 


nisi 


(in  the  case  of  a  constant  outside  pressure,  or  F  —  F(x) 
Eqs.  (8)  will  be  integrated  in  squares). 

]  '  The  values  of  A  and  B  in  the  Asymptotic  .Expreselo: 

|(11)  and  (is)  can  be  determined  by  smsm&st  numerical  cal 


cuiation.  In  case  a 


0.44,  B 


v-jDi  co«=  x ,y;  m 


case  3 >  A  0.166,  B  ^  1.464,  and  cosesI.6. 

The  concentration  of  the  energy  in  the  center  can 
be  investigated  through  the  behavior  of  the  solution  at 
x  Q.  The  energy  per  unit  volume  consists  of  kinetic 
energy,  dissipation  energy  at  the  wavefront,  and  dissipat 
ion  energy  due  to  the  work  of  the  forces  of  plastic  de~ 


these  integrals  at  x  -*♦0  and  s— *■(),  that  the  final  con¬ 
centration  of  energy  in  the  center  takes  place  only  in  th 
case  of  the  collapse  of  a  spherical  layer  of  ineompressib 
Quid  (the  second  integral  in  (14)  is  absent)  under  the 
sondition  nssmk*0  (k  is  arbitrary).  In  other  cases  all 
integrals  In  (14)  converge,  and  this  means  a  decrease  in 
the  stored  energy  at  the  center. 


In  analogy  to  the  above,  the  problem  can  be  solved 
bv  assuming  a  variable  compression  of  matter,  depending 
Up  the  amplitude  of  the  pressure  at  the-  wave.*  r  on  «//*-/* 
Wowever,  ¥a  'in  /3,  V  ^  would  be  necessary  here  to  soivo 
(an  integral  differential  equation  for  tne  vexocuW  of  - 
front  of  the  shock  >»a.ve « 
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